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of zero-range processes and related models. After revisiting some conceptual aspects 
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Introduction 

The models discussed in these notes are simplified models of physical reality. Yet, besides 
the fact that they play an important role in the elucidation of conceptual problems of 
statistical mechanics and probability theory, they are instrumental in the understanding 
of a variety of complex physical situations. 

The aim of these notes is a description of the statics and dynamics of zero-range 
processes (ZRP) [1] and of related models. We present a review of the subject, coming 
back on some of its conceptual aspects. We restrict all discussions to homogeneous 
models where all sites are equivalent. Before commencing, we summarise in a few words 
the main organisation of the text. 

In Part I: Statics (sections 1-6), we first show that ZRP are special members of a 
class of stochastic processes which have the property that their stationary measures are 
known and have a product structure. The probability of a configuration of the system 
is given by the Boltzmann formula for an equilibrium urn model with independent sites. 
Reversibility (for symmetric dynamics) and pairwise balance (for asymmetric dynamics) 
are inherently related to the structure of the stationary measure. Generalisations to 
multiple-species ZRP are then addressed. The properties of the stationary measure of 
ZRP leading to a phase transition between a fiuid phase and a condensed phase are 
finally briefiy reviewed, as a preparation for the second part of these notes. 

The stochastic nature of ZRP is fully revealed by the study of their dynamics. This 
is the subject of Part II (sections 7-9). We first address the nonstationary dynamical 
behaviour of the system when it evolves from a random initial disordered configuration 
to its stationary state. Then we investigate some aspects of its stationary dynamics, 
when the system fiuctuates in its stationary state. In both cases the model used is that 
giving rise to condensation. 

Part I: Statics (sections 1-6) 
1. Dynamical urn models and zero-range processes 

1.1. Dynamical urn models 

We name dynamical urn model (DUM) the following stochastic process. Consider a 
finite connected graph, made of M sites (or urns), on which A^ particles are distributed. 
The occupation Ni{t) of site i [i = 1, . . . , M) is a random variable, and the total number 
of particles 

M 
1=1 

is conserved in time. The model is defined by dynamical rules describing how particles 
hop from site to site. An elementary step of the dynamics consists in choosing a 
departure site d and an arrival site a connected to site d, and in transferring one of 



From urn models to zero-range processes: statics and dynamics 3 

the particles present on site d to site a. This process takes place with rate Wk^i per unit 
time, depending on the occupations both of the departure site, k = N^ ^ 0, and of the 
arrival site, / = Nat- 

On the complete graph (i.e., in the mean-field geometry), all sites are connected, 
i.e., sites d and a are chosen independently at random. On finite-dimensional lattices, 
site a is chosen among the first neighbours of site d. In one dimension, site a is chosen to 
be the right neighbour of site d with probability p, or its left neighbour with probability 
q = 1 — p. In the following we consider the one- dimensional symmetric dynamics, 
corresponding to p = 1/2, and the general asymmetric one, corresponding to p ^ 1/2, 
both with periodic boundary conditions. 

A configuration of the system is specified by the occupation numbers Ni{t), i.e., a 
complete knowledge of its dynamics involves the determination of V{Ni^N2, . . . , Nm), 
the probability of finding the system in a given configuration at time t. 

The process can be pictorially viewed in terms of colonies and migration. The 
sites are the colonies, or cities. An individual leaves its colony for another one, with 
a rate Wk^i which depends on the number of members present in both the departure 
and the arrival colonies. Thus, for example, the philanthrope is characterized by a rate 
decreasing with k and increasing with /, the misanthrope by the converse. 

1.2. Zero-range processes 

Zero-range processes are just particular cases of DUM, with the additional restriction 
that the rate Wk^i only depends on the occupation of the departure site: 

Wk,i = Uk. 

This simple restriction is enough to lead to a remarkable property of the stationary 
probability [1, 2]. Indeed, the probability of a configuration of the system is equal to 

2 M 

P(iVi,...,Ar^,) = - nm, (1-1) 

where it is understood that J2Ni = N, and where the factor pk = PNi=k satisfies the 
relation 

UkPk=Pk-i, (1-2) 

which leads to the explicit form 

1 , , 

Po = 1, Pk = • (1-3) 

Ui. . .Uk 

The normalisation factor, hereafter refered to as the partition function, reads 

Zm,n = E ■ ■ ■ E P^i ■ ■ -P^m 4E Ni, n) . (1.4) 

Ni Nm \ i / 

J Throughout this text wc use the notation Ni for the random occupation of site i, and k (an integer) 
for the value taken by this random variable. 
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One important observation to make is that the stationary measure is insensitive to the 
bias. 

These results can be proved by inspection. The master equation at stationarity 
reads 

0=Y. M{C\C')V{C') - J2 M{C"\C)V{C), (1.5) 

C'^C C'j^C 

where C = {Ni, . . . , Nm}, and M{C\C') is the transition rate from C to C. Consider a 
system of M = 3 sites for simphcity. At stationarity the master equation reads exphcitly 

p [V{Ni + 1, iV2 - 1, Ns) UN,+i{l - 5{N2, 0)) + c.p. + c.p.] 
+ q [P(iVi + 1, iV2, iVs - 1) UN,+iil - 5(iV3, 0)) + c.p. + c.p.] 
= P(iVi,iV2,iV3)[n^^(l-5(iVi,0)) + c.p. + c.p.], (1.6) 

where c.p. stands for circular permutation. Carrying the product form (1.1) into the 
equation, and using (1.2), satisfies the master equation. This is the unique solution of 
the problem. 

The cancelation of terms in the equation occurs by pair. Pairs correspond to terms 
bearing the same p (respectively q) factor, and the same 1 — b{Ni^ 0) factor. Hence we 
have for example 

pV{N^ + 1, iV2 - 1, iVg) M;vi+i = pP(iVi, iV2, iVs) Wiv„ (1.7) 

i.e., using the product form (1.1), with Ni = k and N2 = I, 

Pk+iPi-iUk+i=PkPiUi, (1.8) 

which is precisely the relation that leads to (1.2). 

It is interesting to emphasize the interpretation of (1.7), or (1.8). Consider the 
following configurations: 

C ={N,,N2,N,), 

C =(iVi + l,iV2-l,iV3), 

C"=(iVi,iV2-l,iV3 + l), 

and the corresponding rates 

M{C\C')=puN,+i, 

M{C'\C)=quN„ 

M{C"\C) =puN2- 

In general, i.e. for a general value of < p < 1, (1.7) reads 

M{C\C')V{C') = M{C"\C)V{C). 

This is a condition for pairwise balance [3]. It expresses the equality between the 
probability fluxes flowing from C to C, and from C to C". In the particular case where 
p = 1/2, or more generally when the dynamics is symmetric, then (1.7) becomes the 
condition for detailed balance 

M{C\C')V{C') = M{C'\C)ViC). 
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1.3. Equilibrium urn models with independent sites 

We now adopt a completely different point of view. We consider equilibrium urn 
models with independent sites, on which a dynamics is then defined, in such a way 
that equilibrium is recovered at long times. 

As above, we consider a finite connected graph, made of M sites (or urns), on which 
N particles are distributed. The number of particles on site i is the random variable 
Ni, with Yl^Ni = N. The total energy of the system is defined as the sum 

M 

E(iVi,...,iVM) = E^TO- 

Let 

p^^ = e-^^m (1.9) 

be the unnormalized Boltzmann weight attached to site i. Then, clearly, the probability 
of a configuration of the system is given by the product form (1.1), and Zm,n appears 
as the usual partition function for this statistical mechanical system. 

We now define a dynamics for this model, such that equilibrium is attained in the 
limit of long times. We therefore choose a rule obeying detailed balance for the move 
of a particle. This implies that the dynamics should be symmetric. Restricting to the 
one-dimensional case, [p = 1/2), ii Ni = k and A^i±i = I, we have 

PkPi Wk,i = Pk^iPi+i Wi+i^k~i, (1.10) 

which expresses the probability balance between the configurations {A^^ = k,Na = 1} 
where {d = i,a = i ± 1), and {N^ = I + l,Na = k — 1} where {d = i ± l,a = i). It 
applies as well to the case of the complete graph. For example, with the Metropolis 
rule, the move is allowed with probability min(l,exp(— /^Ai?)), where AE is the change 
in energy due to the move. 

Let us mention two well-studied models in this class: the backgammon model 
and the zeta-urn model, that we briefly describe. The backgammon model is a simple 
example of a system which exhibits slow relaxation due to entropy barriers [4, 5]. The 
following choice of an energy function is done: 

E{N,) = -6{N,,0). 

The statics of this model is trivial. Its interest lies in its dynamical behaviour. The 
dynamics of the model has been thoroughly studied in the mean-field geometry, with 
Metropolis dynamics, and with the additional rule that a particle (instead of a site) is 
chosen at random. The rate for the Metropolis rule reads 

u/ • (^ Pk~iPi+i \ 
Wk,i = mm 1, . 

V PkPi J 
From (1.9), we have 

Po = e^, pk = 1, {k> 1), 
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and therefore H4,o = 6~^ for any k > 1, and W^^i = 1 otherwise, or in compact form: 
Wk,i = 1 + (e-^ - l)6i,oil - Sk,i) {k > 0). 

As can be read on this expression, at fow temperature increasing the number of empty 
sites is not favoured. The total energy is indeed equal to minus the number of empty 
sites, so that particles tend to condensate in fewer and fewer sites as times passes, at 
least at low temperature. 

The static zeta urn model has energy function 

E{N,) = \n{Ni + 1), (1.11) 

hence 

1 
Pk 



The model was initially introduced as a mean-field model of discretized quantum 
gravity [6]. Its dynamics was subsequently defined and investigated in the mean-field 
geometry with heat-bath dynamics [7, 8]. 

If instead, the transfer rate is taken to be that of a ZRP, with Wk^i = Uk, where 
Uk = Pk-i/Pk, the universal properties of the dynamics of the zeta urn model are not 
changed [9]. We thus get 

/ 1\^ (3 

The model is therefore in the same universality class as the ZRP with condensation 
studied in the rest of this text, and defined with the rate Uk = I + b/k. The parameter 
b for this model can therefore be identified with the inverse temperature. 

To summarise at this point, we have so far encountered two classes of dynamical 
urn models with stationary product measures. On the one hand, ZRP are defined for 
any value of the drive, and are such that the transfer rate Wk^i only depends on k. On 
the other hand, equilibrium urn models with independent sites are defined from the 
start without drive, but the transfer rate has the full dependence in both k and /. A 
natural question to ask is whether there exist models possessing both features, namely 
models with stationary product measure, even when submitted to a drive, and with 
transfer rate Wk,i not restricted to depend only on k. 

1.4- Dynamical urn models with stationary product measure 

We address the question just posed. Given a DUM, what choice of rate Wk,i is compatible 
with a stationary measure of the form (1.1), even if the dynamics is not symmetricl 
Let us restrict to the case of the one-dimensional geometry with asymmetric hops. The 
results are as follows: 

• In the general case, < p < 1, two conditions are imposed on the rate W^^i- The 
first condition is 

PkPi Wk,i = Pk~iPi+i Wi+i^k~i- (1-12) 
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The second condition reads 

Wk,i-Wi,k = Wk,o-Wi^o. (1.13) 

Eq. (1.12) expresses the condition of pairwise balance. 

• In the symmetric case {p = 1/2), the only condition imposed on the transfer rate 
is (1.12), or equivalently (1.10). It expresses the condition of detailed balance. In 
other words, if the stationary measure is a product, it is necessarily an equilibrium 
measure and we are taken back to the situation of section 1.3. 

Let us give the proof. By hypothesis, the stationary probability is given and has 
the product form (1.1), with given p^^. We rewrite the master equation (1.5) as an 
equality between gain and loss terms, after dividing both hand sides by V{C), 

pGr + qGl = pLr + qLi, 

with right and left contributions 



i 

Gr= Yl Wn,+j_+i,n,-i 



Pn,+i+iPn,-i 



We now specialize to the configuration where all sites are empty except for sites i 
and i + 1: 



C = {iVi = 0, . . . , iV,_i = 0, iV, = fc, iV,+i = /, Ni+2 = 0, . . . , iVM = 0}. 



We obtain 



..r PlPk^l , ... Pk+lPl-l\ 

P Wi^k^i h Wk+i.i^i 

V PoPk PkPl ) 

, (..r PlPl-1 , TI^ Pk~lPl + l \ 

V PoPi PkPi ) 

= p{Wk,i + Wi^^) + q{Wi^k + Wk,o)- (1.14) 

Taking /c = 0, (1.14) reduces to 

PiPi~i Wi^i^i = poPi Wifl, (1.15) 

which expresses the probability balance between the configurations {A^^^ = 1, A''^ = / — 1} 
and {Nd = l,Na = 0}. This equality is then used in (1.14) to yield the fundamental 
equation 

p{Wk,i - Wk,o) + q{Wi,k - Wi,o) = 

P (Wk^u-i ^-^^^^ - W,^ + q iw.^.^k-. ^-^^^^ - Wk^ . (1.16) 

V PkPi J \ PkPi J 

From this equation, the two conditions (1.12) and (1.13) are obtained, as shown in the 
appendix. The conditions thus found are necessary. They are also sufficient as one can 
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convince oneself by redoing the reasoning for a generic configuration. The analysis done 
here applies as well to the complete graph, for which the dynamics is symmetric. 

Coming back to the case of a ZRP, condition (1.13) is trivially satisfied, while the 
pairwise balance condition (1.12) yields (1.8), rewritten here for convenience, 

PkPlUk = Pk-lPl+lUl+l- 

The ZRP appears as the minimal model of the class of DUM leading to a product 
measure in the stationary state independent of the asymmetry. It is important to 
realize that this measure is that of an equilibrium urn model with independent sites 
(see section 1.3) and therefore any result on the statics of a ZRP pertains to the field 
of equilibrium statistical mechanics. 

The original work on the question posed in the present section is due to [10]. The 
dynamical urn model described in the present notes is named a misanthrope process 
in [10] because the rates Wk,i considered in this reference are increasing functions of k. 
Yet another presentation, restricted to the ID totally asymmetric case {p = 1) can be 
found in [11]. 

2. A counterexample 

Let us now examine the case where the transfer rate only depends on the occupation of 
the arrival site, 

Wk,i = vi{l-6k,o). (2.1) 

If the dynamics is symmetric, the only constraint to take into account in order to have 
product probability in the stationary state is the detailed balance condition (1.12), which 
reads here pk+iPi^i fz-i = PkPi Vk- The relation pi_i vi^i = pi follows, which determines 
the measure fully. However, if the dynamics is not symmetric, (1.13) is violated by (2.1), 
which rules out the possibility of stationary product measure for this case. 

Let us illustrate the difficulty on the simple case of a system of M = 3 sites. The 
stationary master equation reads 

p [P(iVi + 1, iV2 - 1, iVa) VN,^^{1 - 5{N2, 0)) + c.p. + c.p.] 
+ q [P(iVi + 1, iV2, iVg - 1) i;jv3_i(l - <5(iV3, 0)) + c.p. + c.p.] 
= pP(iVi,iV2,iV3)[^^^i(l-5(iV3,0)) + c.p. + c.p.] 
+ gP(iVi,iV2,iV3)[t'^i(l -'^(iV2,0)) + c.p. + c.p.]. 

There is no way of pairing the terms in the master equation to obtain their mutual 
cancellation iip ^ 1/2, while this is possible for p = 1/2. More generally, the stationary 
probability is unknown for the asymmetric process (for arbitrary system size M) [12]. 

3. Two-species ZRP: conditions for product measure 

A simple generalization of the ZRP defined so far consists in considering two (or more 
generally n) coexisting species on each site [13, 14], named particles of type A and B 
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respectively. The hopping rates for A and B particles only depend on the occupations 
of the departure site: Nf- = k, Nj^ = I. They are respectively denoted by u^^i and Vk^i- 
The new fact is that the condition for product stationary measure imposes a constraint 
on the rates Uk^i and Vk,i [13, 14], given by equation (3.4). 

We revisit this problem, keeping the line of thought followed for the (single-species) 
ZRP in section 1. We want to show that, as was the case for the single-species ZRP, for 
the two-species ZRP satisfying (3.4) the following properties come together: 

• The stationary probability is a product and is insensitive to the bias. It is the 
stationary probability of an equilibrium urn model with independent sites. 

• If the dynamics is symmetric, the process is reversible, i.e. satisfies detailed balance, 
otherwise, in the presence of a bias, pairwise balance holds. 

3.1. Equilibrium urn models with independent sites 

Let us first consider an equilibrium urn model for two species with independent sites. 
A configuration of the system is denoted by C = {Ni, . . . , Nm}, where Ni = {Nf-, N^). 
The energy is given by the sum 

E(c) = E^(M). 

i 

The Boltzmann weight reads 

nC) = ^ UPn.^ (3-1) 

where A^"^ and A^^ are respectively the total number of A and B particles, and Z the 
partition function. A dynamics yielding this equilibrium measure should fulfill detailed 
balance. We restrict the rates to depend only on the departure site. With for the 
departure site: N^ = k, N^ = I, and the arrival site: N^ = m, N^ = n, we must 
impose 

Pk,lPm,nUkJ = Pk~l,lPm+l,nUm+l,n 

Pk,lPm,nVk,l = Pk,l-lPm,n+lVm,n+l, (3.2) 

hence, 

Uk,iPk,i = Pk-i,h Vk,iPk,i= Pk,i-i- (3.3) 

These relations generalise (1.2). Consideration of the two possible paths leading from 
Pk,i to pk-i,i-i, using (3.4), imposes a "gauge" condition on the rates: 

Uk,iVk-i,i = Vk,iUk,i-i. (3.4) 

3.2. Product measure 

We can now proceed as for the single-species ZRP. We claim that, even in the presence 
of a bias, (3.1), (3.3), hold. The proof is by inspection: (3.1) carried into the master 
equation of the process is seen to be a solution if (3.2), or (3.3), hold. The constraint 
(3.4) follows. See [13, 14]. 
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3. 3. Reversibility implies stationary product measure 

Finally we show by a direct route that (3.4) is a consequence of reversibility, when the 
dynamics is symmetric. We use the Kolmogorov condition, a necessary and sufficient 
condition for the reversibility of a Markov process (i.e., for detailed balance to hold), 
which states that the product of rates along any cycle in the state space of the process 
and for the reverse cycle should be equal [15, 16]. Consider the configuration 

(iVi = (A;,0,iV2 = (m,n),...). 

We consider the following cycle in the space of states of the process 

(/c, /; m, ra) — > {k,l — l;m,n + 1) -^ {k + 1,1 — l;m — l,n + 1) 
— s> (fc + 1, /; m — 1, n) — > {k, I; m, n). 

For the cycle considered above, the Kolmogorov condition yields 

'Vk,l "^fc+l.i _ Um,n Vm,n+1 
Vk+l,lUk+l,l~l Um,n+1 Vm-l,n+l 

This condition is satisfied if and only if (3.4) holds. 

3.4- An example of a two-species ZRP with non product stationary measure 
Consider the ZRP defined by the following rates [17] 

Uk,l = 1 + -, Vk,l = 1 + -. 

These rates violate (3.4), and therefore, as explained above, this process violates time 
reversal symmetry even in the absence of a bias. The study of the stationary properties 
of the model is addressed in [18]. 

4. Two extreme cases 

4.I. The case of two sites 

We come back to the case of a general dynamical urn model, with one species, where 
now the number of sites is M = 2. This case is interesting for several reasons. Firstly 
the model stands by itself, for instance the Ehrenfest urn model belongs to this class, 
as shown below. Secondly, it illustrates some aspects of the general theory for a system 
of arbitrary size M. Finally, it relates to the other case considered in this section, a 
thermodynamic system on the complete graph, the master equation of which is formally 
that of a two-site system. 

Since N2 = N — Ni, a configuration of the system is entirely defined by the 
occupation of site 1, A'^i, and the hopping rate only depends on one variable: Wk^i = Uk- 
Let us denote the occupation probability of site 1, i.e., the probability of a configuration 
of the system, by 

fk{t) = V{N,{t) = k). 
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It obeys the master equation 
dfkit) 



dt 

dm 

dt 
dfNit) 



Hk+i fk+i + Afc-i fk-i - [i^k + Afc)/fc (1 < A; < iV - 1), 
/ii/i-Ao/o, (4.1) 

Aat-i Jn-i — (^nIn, 



dt 

where A^ and Hk are respectively the rate at which a particle enters site 1, coming from 
site 2, or leaves site 1 for site 2: 

Afc = UN-k, f^k = Uk- 

The equations for k = or k = N are special, since uq = 0. The above equations 
describe a biased random walk on the interval (0,iV), with reflecting boundaries at 
and N, the position of the walker being the random variable Ni{t), i.e., the number of 
particles on site 1. 

The time-independent solution to (4.1) satisfies 

/^fc+l /fc+l,eq ~ Afc fk,eq = . . . = /il /l,cq ~ Aq /o,cq = 0, 

which yields the detailed balance condition at equilibrium 

fJ'k+l fk+l,oq = Afc /fc,eq. (4.2) 

From this equation it is easy to obtain 

N 

r PkPN-k ry _ ^^ 

jA:,eq — — 7} ) ^2,Ar — / ,VkVN-k-, 

^2,7V fc=o 

where the pfc are given by (1.3). These expressions are special instances of eqs. (1.1) and 
(1.4) which hold for the general case. Elements on the dynamics of the two-site model 
can be found in [19]. 

Remark. This process is equivalent to the historical Ehrenfest model [20, 21], defined 
as follows. Consider N particles, labeled from 1 to A^, which are distributed in two urns 
(sites) . At random times, given by a Poisson process with unit rate, a particle is chosen 
at random (i.e., an integer between 1 and A^ is chosen at random), and moved from the 
site on which it is to the other site. The master equation reads 

'-f^-'^M).^^^f.-dt)-Mt). (4.3) 

Indeed, a move of a particle from site number 1 to site number 2 (resp. from site number 2 
to site number 1) occurs with a rate k/N (resp. (A^ — k)/N) per unit time. 

Note that the rule of choosing a labeled particle is different from the rule adopted 
above for dynamical urn models (as was already the case for the backgammon model). 
Yet we can describe this model as a 2-site dynamical urn model, by taking Uk = k 
(dropping the factor A^ which enters the scale of time). Then pk = ^/k\, and the 
distribution of particles amongst the two sites is binomial, 

fk,eq = 2-''(j] ik = 0,...,N), (4.4) 

as is well known for the Ehrenfest model. 
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4-2. A thermodynamic system on the complete graph 

In the mean-field geometry, for a tliermodynamic system, tlie temporal evolution of the 
occupation probability fk(t) is given by the master equation 

dfk{t) 



dt 

dfojt) 

dt 



/ifc+i fk+i + Afe_i /fc-i - (/ifc + Afc)/fc {k > 1), 

/Ui /i - Ao/o, (4.5) 



where 



Hk = Uk, ik>0), \k = ^uifi=ut, {k>0). (4.6) 

1=1 

These are respectively the rates at which a particle leaves site 1, or arrives on this site. 

In other words, on the complete graph, all sites other than site 1 play the role of a single 

site from which particles are emitted with rate Ut, and therefore (4.5) is formally similar 

to the master equation (4.1) for a system of two sites. In the present case this set of 

equations is non linear because Ut is itself a function of the fk{t). 

In the stationary state the detailed balance condition (4.2) reads 

Jk+l,eq _ Afc _ Meq 



yielding 



/fc,eq /^fc+1 Uk+1 



_ Ao . . . Afc_i 

Jk,eq JO,cqi 

/il.../ifc 



where /o,eq is fixed by normalisation. Hence 



l,eq 



Jk,eq — ^oo „ -fc ' ('^- ') 

eq 



PkU. 

with the pk given by (1.3). This expression is a particular instance of the general case 
(5.9). 



5. Statics of ZRP: Fundamental properties 

We collect here the results found so far concerning single-species ZRP's. A ZRP is a 
dynamical urn model, for which the rate of transfer of a particle, Uk, only depends on 
the occupation of the departure site, k. The stationary state of a ZRP is that of an 
equilibrium urn model with independent sites: the probability of a configuration of the 
system is (independently of the asymmetry) 

^ M 

p(iVi,...,iVM) = ^ — Hm, (5.1) 

with partition function 

Zm,n = E ■ ■ • E P^i ■ ■ -P^m 4E N^, n] . (5.2) 
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The factor p^v. obeys the pairwise balance condition (1.12), i.e., PkPiUk = Pk-iPi+iUi+i, 
and hence 

PkUk = Pk-l, 

which gives the exphcit form of pk (for Uk given) 

1 , , 

Po = 1, Pk = • (5.3) 

Ui. ..Uk 

The value given to po is arbitrary. The energy function associated to the underlying 
equilibrium urn model mentioned above is defined using eq. (1.9). 
The partition function Zm,n obeys the recursion formula 

TV 

Zm,N = 2^PkZM-l,N-k- (5.4) 

fe=0 

This ensures that the stationary single-site occupation probability 

f TllAT 1\ PkZM-l,N-k /p. p.x 

/fc,st = V{Ni = k) = (5.5) 

is normalised. We have 

N 

Zo,N = Sn,0, Zi^N = Pn, Z^^N = / ,PkPN-k, (5.6) 

fc=0 

and so on. Using an integral representation of the Kronecker delta function, 

6(m,n) = <p ttZ^-, 

we obtain 

d^ 



^M,N 



-P{zY\ (5.7) 



l-KXZ^ 

where the generating series of the weights pk reads 
P{z) = Y.Pkz\ 

k>0 

In other words, Zm,n is the coefficient of z'^ in P{zY'^ . Static properties of the ZRP are 
therefore entirely encoded in this series. 

In the thermodynamic limit (M — > oo at fixed density N/M = p), the free energy 
per site, 

J^ = - lim -— In Zm,n, 

can be obtained by evaluating the contour integral in (5.7) by the saddle-point method. 
The saddle-point value zq depends on the density p through the equation 

zoP'{zo) , , 

-n^ = "■ ''^" 

The free energy per site is JF = pln^o ~ fi^-P(2^o); and the stationary occupation 
probability reads 

k 

^'^'^ = IS) • ^^-^^ 
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Eq. (5.8) can be rewritten as 

{N^) = Y,kfk,st = p. (5.10) 

k 

Note that the function 

is increasing with zq because 

dp{zo) ,. ,. 
Zo = vai Ni. 

dzo 
Finally the stationary average rate reads 

- fA/T AT\ I \ V^ t V^ PkZM-l,N-k Zm,N-1 /r-ll\ 

Ust{M,N) = {un^} = 2^Ukfk,st = l^Uk = . (5.11) 

k k ^AI,N ^M,N 

In the thermodynamic limit, we have Ugt = Zq (defined in (5.8) above). The expression 
(4.7) found for the case of the complete graph in the thermodynamic limit is a particular 
example of (5.9). 

6. Statics of ZRP: Examples and the phenomenon of condensation 

We illustrate through examples the considerations of the previous section. In particular 
we discuss the possible solutions of eq. (5.8) (or (5.10)). Two possible situations can 
arise. Either p{zq) is allowed to increase without bounds, in which case the equation 
has a solution in zq for any value of p. Or p{zq) reaches a maximal value, Pc, in which 
case the equation has no solution if p > pc. 

6.1. Two simple examples 

Let Uk = k. This model can be seen as a multi-urn generalisation of the Ehrenfest 
model. We have P{z) = e^. The radius of convergence of this series is infinite. Hence 
eq. (5.8) has a solution for any value of p: p(zo) = ^o, hence zq = p, and 

h t = e-"^ 

J k,st " , • ; 

which is a Poisson distribution. The fast decay of the distribution is characteristic of 

an homogeneous fiuid phase. 

As a second example let m^ = 1. Then P{z) = 1/(1 — z). The partition function of 

a finite system is 

_/M + N-l\ 
Zm,n -iy ^ j. 

The radius of convergence of P{z) is equal to 1. At this maximal allowed value of 
z, p{z) = z/{l — z) is infinite. Therefore (5.8) has a solution for any value of p: 
Zq = p/{l + p), and finally 

f - ^ I P \ 

J fc,st 



1 + p \l + pj 

The system is again in a fiuid phase. 
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6.2. The canonical example for the phenomenon of condensation 

We consider the ZRP with transfer rate 

b 

This case, and closely related models, have been studied in various references [6, 22, 
23, 24, 9, 25, 7, 8, 26, 27]. We follow here the approach and notations of [19]. For this 
choice of rate, 

Pk ... 



r{k + b+l) Jo ^ ' k^ ' 

/•I h(l - u)''~^ 

P{z) = / du \ ^— = 2Fi(l, 1; 6 + 1; ^), (6.1) 

JO I — zu 

where 2F1 is the hypergeometric function. The function P{z) has a branch cut at z = 1, 

with a singular part of the forni§ 

P,,{z)^AP{l){l-zf-\ 

so that P{z) is only differentiable n = Int(6) — 1 many times at 2; = 1: 

P{Z) ^ P(l) + (1 - ^) P'(l) + • ■ ■ + ^ p^P(")(l) + Psg(^). 

The following values are of interest: 



0—1 Sm TTO 

b Ab 



S.2) 



(6-l)(6-2)' '' {b - l){b - 2){b - 3)' 

For b <2, p(l) is infinite. The system is in a fluid phase: 

/fc,st~A;-''e-'=|i°^»l (6.3) 

For 6 > 2, p(l) is finite. The system has a continuous phase transition at a finite 
critical density 

P'jl) _ Ekkpk 1 

^' P(l) E.P. &-2' 
such that the saddle point zq reaches the singular point z = 1. This critical density 
separates a fluid phase (p < pc) and a condensed phase (p > pc). 

Fluid phase (p < pc). The equation (5.8) has a solution for any p < pc- The single site 
probability has the form (6.3). 

§ Whenever 6 = n > 2 is an integer, the amphtudc A diverges. The singular part of the generating 
series is of the form Psg(^) ~ n(— 1)"(1 — z)"~^ ln(l — z). 
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Critical density (p = pc) ■ The occupation probability 

p, ^ Ab-i)m 

Jk,st — p^x ~ ^j lu.-t; 

falls off as a power-law in the thermodynamic limit. The critical free energy reads 

b 



Tr = -InPd) = -In 



6-1 
The second moment of the occupation probability, 

is convergent for 6 > 3 (regime of normal fluctuations), and divergent for 2 < 6 < 3 
(regime of anomalous fluctuations). 

Condensed phase (p > p^). A large and finite system in the condensed phase essentially 
consists of a uniform critical background, containing on average N^. = Mp^ particles, 
and of a macroscopic condensate, containing on average A = N — N^ = M{p — pc) excess 
particles with respect to the critical state. 

The occupation probability fk^st accordingly splits into two main contributions [24]. 
The contribution of the critical background, corresponding to small values of the 
occupation [k <^ M), is approximately given by (6.4). The contribution of the 
condensate shows up as a hump located around k = A. The hump is a Gaussian 
whose width scales as M^/^ whenever pc is finite, i.e., for b > 3, whereas it has power- 
law tails and a larger width, scaling as M^^^''~^\ in the regime of anomalous fluctuations 
(2 < 6 < 3). The weight of the condensate probability hump is approximately 1/M, 
in accord with the picture that the system typically contains a well-defined condensate 
located on a single site at any given time. 

6.3. Rate Uk = 1 + a/k'^ : Stretched- exponential critical behaviour 

Consider the ZRP with hopping rate [22, 26, 19] 

Uk = l + —, (6.6) 

where a is an arbitrary exponent. The situation of interest corresponds to < cr < 1. 
Equation (1.3) leads to the estimate 

Pfc ~exp f-a^ — J ~exp(^-^— -A;^-'"j . (6.7) 

The generating series P{z) has an essential singularity at 2; = 1 with an exponentially 
small discontinuity. The critical density 
_P'il) _ T.kkpk 

^^ P{l) EkPk 

is finite. The occupation probability at the critical density, fk^st = Pkl P{^)i decays as a 
stretched exponential law. 



From urn models to zero-range processes: statics and dynamics 17 

Part II: Dynamics (sections 7-9) 

7. Zero-range processes: nonstationary dynamics (I) 

The question is to determine the temporal evolution of the system starting from a 
random disordered initial condition. Here we study the dynamics of the class of ZRP 
giving rise to a condensation transition in their stationary state. For simplicity we will 
choose the hopping rate 
, b 

We address the question first in the fully connected geometry. 

The same question can be addressed for dynamical urn models (see e.g. [5]). The 
analysis that follows [9] , as well as that contained in the next section, are essentially the 
same as that performed for the zeta-urn model [7, 8]. 

7.1. Dynamics on the complete graph 

We wish to determine the temporal evolution of the occupation probabihty /^(t). 
Conservation of probability and of density yields 

oo 

E/^W =1' (7-1) 

fc=0 

oo 

T.kfk{t)=p, (7.2) 

fc=i 

where we have taken the thermodynamic limit A^ -^ oo, M -^ oo, with fixed density 

p = N/M. We consider a system with Poissonian initial distribution of occupation 

probabilities, 

MO) = e-|^, 

i.e., such that initially particles are distributed at random amongst sites. 

Since the equations (4.5) are non linear they have no explicit solution in closed 
form. Yet one can extract from them an analytical description of the dynamics of the 
system at long times, both in the condensed phase, and at criticality The structure of 
the reasoning borrows to former studies on urn models [7, 8]. (For a review, see [5].) 

As we show below, there exists two different regimes in the evolution of the system, 
both in the condensed phase or at criticality, which we study successively. 

(a) N onequilihrium dynamics of condensation (p > pc) 
Since -Ugq = 1, we set, for large times, 

Ut^l + Ast, (7.3) 

where the small time scale et is to be determined, and A is an unknown amplitude. 
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Regime I: k fixed, t large. For t large enough, sites empty {uk) faster than they fill (mj). 
In this regime there is convergence to equilibrium, hence we set 

fk{t)^hW^+vuet), (7.4) 

with /fceq given by (6.4), and where the Vk are unknown. This expression carried into 
(4.5) yields the stationary equation /^ = 0, because the derivative /fc, proportional to 
St, is negligible compared to the right-hand side. We thus obtain an equation similar to 
the detailed balance condition: 

/fc+l,eq 1 + Vk+1 £t _ 1 + Aet 

fk,eq '^ + VkEt Uk+l 

Using (6.4) and (1.3), we obtain, at leading order in et, u^+i — v^ = A, and finally 

Vk = Vq + kA. (7.5) 

At this stage, vq and the amplitude A are still to be determined. 

Regime II: k and t are simultaneously large. This is the scaling regime, with scaling 
variable x = kst- Following the treatment of [7, 8], we look for a similarity solution of 
(4.5) of the form 

h{t)^{p-p.)eU{x). (7.6) 

We thus obtain for g{x) the linear differential equation 

with St ~ t~2. This is precisely the differential equation found in [7, 8], for the zeta- 
urn model. The amplitude A can be determined by the fact that the equation has an 
acceptable solution g{x) vanishing as x ^ and a; ^ oo [7]. The amplitude A and the 
scaling function g{x) are universal quantities, only depending on the value of h. The 
sum rules (7.1) and (7.2) yield respectively 

/ (lxg{x) = — -, 

Jo Pc- P 

POO 

/ dxxg{x) = 1. 
Jo 

The differential equation above has no closed form solution. However further information 

on the form of the solution g{x) can be found in [7, 8]. 

An intuitive description of the dynamics of condensation in the scaling regime is 

as follows. The typical occupancy A;cond of the sites making the condensate scales as tz . 

The total number of particles in the condensate is equal to M{p — p^), the remaining 

M Pc lying in the fluid. Therefore the number of sites belonging to the condensate scales 

as M{p — pc)t~"^. 
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(b) N on equilibrium critical dynamics (p = Pc) 
The analysis follows closely that done in [8]. We set 
ut^ 1 + Aet, 

with Et = t~'^, where the exponent u is to be determined, and we consider the same 
two regimes as above. In regime I, we still set (7.4) for fk{t). The reasoning leading to 
the relationship v^ = Vq + k A (see (7.5)) is still valid here. In regime II, we look for a 
similarity solution to (4.5) of the form 

fkit) ^ fk,cqgcix) x = kt~^. (7.7) 

Indeed, for any large but finite time t, the system looks critical, i.e., the occupation 
probabilities fk{t) have essentially converged toward their equilibrium values (6.4), for 
k -C t^/^, while for k ^ t^/^ the system still looks disordered. The fk{t) are expected to 
fall off very fast, which is confirmed by the following analysis. 

The sum rules (7.1) and (7.2) lead respectively to the following equations, provided 
that 6 > 3, 





vo + Apc 


= 0, 








(7.8) 




r" {vop. 


., + Ap,) =r(^-2)/2(^_ 


-i)r{b) r duu^~\i- 

Jo 


-9c{u)), 


(7.9) 


where fie - 


= J2k'^ /fc,cq is given in eq. (6.5). Equation (7.9) 


fixes the value of u: 






u = {b- 


-2)/2. 








(7.10) 


The differential equati 
reads 


ion obeyed by gdx) is 


obtained by 


carrying 


(7.7) into 


(4.5). It 


-1- ^^\.AJ\^A-KJ 


9cix) + 


/x b\ 










the solution of which is, with gdO) = 1, 












9cix) = 


^"' rdvv'e-y'/' 








(7.11) 




TCf)L ^^^' • 



The fall-off of gdx) for x ^ 1 is very fast: gdx) ~ exp(— x^/4), hence fk{t) ~ 
exp(— /c^/4t). We finally obtain 

fie- pi Jo 6-1 \2J 

Let us mention that for any hopping rate of the form u^ ^ 1 + b/k, the scaling 
functions, g{x) in the condensed phase (more precisely: g{x)/{p — pc)), and gdx) at 
criticality, are universal. In both cases the scaling variable is x = kt~^^'^. The critical 
density pc, and, as a consequence, any quantity depending on pc, such as the amplitude 
Vq, are non universal, with values depending on the precise definition of Uk- As noted 
above, the amplitude A is a universal quantity in the condensed phase. 
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7.2. Late stages of the dynamics and the case of one dimension 

As mentioned above, in the first stage of tlie dynamics, in tlie MF geometry, tlie number 
of most populated sites decays as M/t^/^. Hence, after a time of order M^, the system 
contains a finite number of highly populated sites, i.e., condensate precursors. 

The late stage of the non-stationary dynamics, where all but one of the precursors 
die out, is thus expected to also last a length of time of the order of the diffusive 
timescale M^. This is substantiated by numerical simulations in [25]. Another argument 
is presented in section 9.4. The whole non-stationary process of the formation of the 
condensate is therefore characterised by a single timescale 

r2 



' non— st 



M\ 



The same results hold for the IDAS case. The analysis relies upon numerical work 
or heuristic and scaling arguments [9, 25]. 

A similar scenario holds in the IDS geometry, the only difference being that Tnon-st 
now scales as M^. The shift of the dynamical exponent by one unit in the IDS geometry 
has a common origin [9, 25]: it stems from the Gambler's ruin problem [28]. An 
analogous phenomenon is encountered for example in the coarsening law for the domain 
growth, and in the motion of a tagged particle, in ID Kawasaki dynamics [29]. 

We refer to the original references for further results (scaling functions, critical case, 
etc.). 

8. Zero-range processes: nonequilibrium dynamics (II) 

So far we considered the dynamics of one-time quantities, related to the random variable 
Ni{t). We now explore another facet of the nonequilibrium dynamics of the ZRP with 
hopping rate Uk = 1 + b/k, namely the two-time nonstationary aspects of its dynamics. 
This essentially means that any function of the two times depends on both times, instead 
of depending on their difference, which would be the case at stationarity. The situation 
here is analogous to that encountered when a ferromagnetic spin system is quenched 
from a high temperature, corresponding to an initial disordered configuration, to a lower 
temperature, T <Tc [30, 31]. 

We consider the same ZRP as in the previous section, on the complete graph, in 
the thermodynamic limit. The system relaxes from a nonequilibrium initial condition 
towards equilibrium. In order to characterize the fluctuations of the local density 
of particles, Ni(t), around its mean {Ni{t)) = p, we study its associated two-time 
correlation and response functions, and fluctuation-dissipation ratio. 

8. 1 . General framework 

The connected two-time correlation function of the density between time s (waiting 
time) and time t (observation time), with s < t, is defined as 

C{t,s) = {N,{s)N,{t))~p\ 
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It can be rewritten as 



C{t,s) = Y.k-fk{t,s)-p\ 

k>l 

where the function 7fc(t, s) is defined by 

7k{t,s) = Y.jf,is)V{N,it) = k I iVi(s) = j} 

with the initial value at t = s 

7fc(s,s) = kfk{s). 

Its temporal evolution for t > s is given by the master equation (4.5): 
d'ykit,s) 



dt 



yUfc+i 7fc+i + -^fc-i 7fc-i - (/Wfc + hhk {k> 1), 

Q /ii7i-Ao7o- (8.1) 

The rates A^ and fik are defined in (4.6). The rate A^ only depends on the fk{t), 
hence (8.1) are linear equations for the •jkit, s). 

The local response function measures the infiuence on the mean density on site 
number 1 of a perturbation in the canonically conjugate variable, i.e., the local chemical 
potential acting on the same site. Suppose that site number 1 is subjected to a small 
time-dependent chemical potential ai(t), so that the total reduced energy of the system 
(see section 5) is now 

M 

PE{{N,}) =Y.PE{N,) + ai{t)Ni. 

i=l 

The mean density on site number 1 reads 

{Ni{t)) =p+ f dsR{t,s)ai{s) + ---, 
Jo 

where only the term linear in a{s) is written explicitly. The kernel of the linear response 

is the two-time response function 

The temporal evolution of this function is given by a master equation similar to (4.5) [8]. 
The zero-range processes that we consider here have a fast convergence towards 
equilibrium, with a finite relaxation time Trdax in their fiuid phase, as is the case for 
a generic statistical-mechanical model in its high-temperature disordered phase. If the 
earlier time exceeds the relaxation time (s ^ Trdax)) the system is at equilibrium. 
One-time quantities take their equilibrium values. Two-time quantities, such as the 
correlation and response functions, are invariant under time translations: 

C{t, s) = Ceq(r), R{t, s) = /2eq(r), (8.2) 

where t = t — s > 0. They are related by the fiuctuation-dissipation theorem 

R.,{r) = -^^. (8.3) 
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In the condensed phase and at criticahty the relaxation time Troiax becomes infinite. 
If the waiting time s and the observation time t are much smaller than Ti-dax; both time- 
translation invariance (8.2) and the fiuctuation-dissipation theorem (8.3) are violated. 
It is convenient [32] to characterize departure from equilibrium by the fiuctuation-dissi- 
pation ratio 

Xit.s)^^^ (8.4) 

ds 
In general, this dimensionless quantity depends on both times s and t and on the 
observable under consideration. It may also exhibit a non-trivial scaling behavior in 
the two-time plane. In all known cases it is observed that 

0<X{t,s) < 1. 

8.2. Application: ZRP with condensation (uk = 1 + h/k) 

N onequilihrium critical dynamics {p = pc) . Let us first note that the variance of the 
population of site number 1 converges to its equilibrium value C^q = fic — pl ^^ ^ power 
law: 

23-6^-(fe-3)/2 

The derivation of the behaviour of the two-time density correlation and response 
functions is the same as in [8]. In the nonequilibrium scaling regime {s,t ^ 1), one 
finds 

C(t,s)^S-(''-3)/2$(x), 

dCit, s) 



s 



(^'i)/2^i(x), (8.6) 



ds 

i?(t,s)^S-(^-l)/2$2(x), 

where 

x = t/s> 1. 

As a consequence, in the scaling regime, the fiuctuation-dissipation ratio X(t, s) only 
depends on x: 

X(t,s)^A'(x) = |44- 
^i{x) 

The dimensionless scaling function X{x) is universal, and it admits a non-trivial limit 

value in the regime where the two time variables s and t are well separated in the scaling 

regime [33]: 

Xqo = lim lim X(t, s) = A'(oo). 

Explicit expressions for the above scaling functions can be derived, using a spectral 

decomposition in Laguerre polynomials [8]. The limit fiuctuation-dissipation ratio thus 

obtained 

^+1 ., 
Xoo = (b > 3 , 
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lies in an unusually high range (4/5 < X^ < 1) for a critical system. Indeed, statistical- 
mechanical models such as ferromagnets are observed to have < Xq^ < 1/2 at 
their critical point. The upper bound X^o = 1/2, corresponding to the mean-field 
situation [33], is also observed in a range of simpler models [32, 34]. 

The above results illustrate general predictions on nonequilibrium critical dyna- 
mics [35, 34, 33, 31]. The exponent of the waiting time s in the first line of (8.6) already 
appears in (8.5). It is related to the anomalous dimension of the observable under 
consideration, and would read {d — 2 + rj)/ Zc for a d-dimensional ferromagnet, where 
Tj is the equilibrium correlation exponent and z^. the dynamical critical exponent. The 
scaling functions $(x), $1,2(2;) are universal up to an overall multiplicative constant, and 
they obey a common power-law fall-off in x~^l'^ . The latter exponent is not related to 
exponents pertaining to usual equilibrium critical dynamics. It reads —\c/zc = Qc—d/zc 
for a ferromagnet, where Ac is the critical autocorrelation exponent [36] and 6c is the 
critical initial-slip exponent [35]. 

Nonequilibrium dynamics of condensation (p > pc). In the scaling regime, two-time 
quantities are found to scale as [8] 



C(t,s)^(p-Pc)s^/'$(x), {x = t/s) 

^^-(P-Pc)s-/^$1(X), 

X{t,s)^X{x)- ^^^""^ 



(8.7) 



$i(x) 

The scaling functions $(x), $1,2(2;) have finite values, both at coinciding times [x = 1) 
and in the limit of large time separations {x = 00). The limit fiuctuation-dissipation 
ratio Xqo = A" (00) depends continuously on b throughout the condensed phase (fe > 2), 
and vanishes only as 

for b large, which corresponds formally to low temperature, while coarsening systems 
are known [37] to have identically Xoo = throughout their low-temperature phase. In 
figure 1 a summary of the values of Xoo is presented. 

This dynamics is different from the usual phase-ordering dynamics [30]. Indeed, 
when a ferromagnet is quenched below its critical temperature, domain growth and 
phase separation take place in a statistically homogeneous way, at least for an infinite 
system. In the present situation, condensation takes place in a very inhomogeneous 
fashion, since fewer and fewer sites are involved in the process. 

8.3. One dimension 

For both the symmetric and asymmetric cases the response can be defined in the same 
fashion as above. There is no analytical tools at our disposal to compute these functions, 
even in the scaling regime. However, for the symmetric case, the fiuctuation-dissipation 
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Figure 1. Plot of the limit fluctuation-dissipation ratio X^o against b. Upper 
curve: critical point {b > 3, p = pc). Lower curve: condensed phase (6 > 2, 
p > Pc)- Thin dashed lines: continuation of the results to lower values of b. 



still holds at equilibrium, while it should be violated in the stationary state of the 
asymmetric case. 



9. Stationary dynamics of the condensate 

9.1. The question posed 

Consider a ferromagnetic system, an Ising spin system for instance. At equilibrium in 
the low temperature phase, the spin symmetry is spontaneously broken. There are two 
possible equilibrium states, one with positive magnetisation, the other one with negative 
magnetisation. However, if one observes a large but finite system, then as time passes, 
the magnetisation keeps changing sign, the system flipping between the two possible 
equilibrium states. Ergodicity is restored for a finite system. The typical time between 
two flips is exponential in L'^~^ ., where L is the linear size of the system, and d the 
dimension of space. 

A similar situation occurs for in the condensed phase of a ZRP. Here the 
spontaneously broken symmetry is translational invariance. For a large but finite system 
in the stationary state, as time passes, the condensate keeps moving across the system. It 
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Figure 2. Dynamics of the condensate (IDS geometry with 6 = 4, M = 40, 
N = 80). Upper paneh instantaneous number of particles A^max(i) on the most 
populated site. Lower panel: location imax(i) of that site. 



spends long lengths of time on a given site, before suddenly disappearing and reappearing 
on another site. The typical value of these lengths of time defines the characteristic time 
T of the dynamics of the condensate. The aim of this section is to analyse the nature of 
this motion and in particular to characterise how r scales with the system size M. 



9.2. Numerical observations 

An intuitive understanding of the phenomenon is easily gained by performing Monte- 
Carlo simulations. These simulations, done in the three geometries: mean-field (MF), 
one-dimensional asymmetric (IDAS) (p = 1), and one- dimensional symmetric (IDS) 
{p = q = 1/2), lead to a common picture. 

The condensate is immobile for rather long lapses of time; it then performs sudden 
random non-local jumps all over the system, at Poissonian times whose characteristic 
scale grows rapidly with the system size M. Figure 2 illustrates this process for the 
IDS case, for a system of size M = 40, with A^ = 80 particles, i.e., p = 2, and b = 4, 
hence Pc = 1/2. The upper panel shows the track of the instantaneous number of 
particles A'max(i^) on the most populated site. The signal for A^max(i) fluctuates around 
A = M{p — Pc) = 60, the mean size of the condensate. The lower panel shows the label 
^max(^) of that site, i.e., the location of the condensate. The non-local character of the 
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motion of the condensate is clearly visible, whereas the longest lapses of time where the 
condensate stays still give a heuristic measure of the characteristic time r. 

We show in what follows that r ~ M^ for the fully connected geometry and the 
directed case, while r ~ M^"*"^ for the symmetric case. Moving the condensate is 



therefore slower than forming the condensate (rnon-st 
see section 7.2). 



M^, Tnon-st ~ M^ respectively. 



9.3. Theoretical analysis 

All the idea relies on a problem of barrier crossing. Defining the potential as Vk = 
— In/fcst; then a dip in the probability fk^st corresponds to a barrier in the potential. 
The flipping time r is the time to cross the barrier, or the first-passage time from right 
to left. Let us explain these ideas in more detail. 




Figure 3. Logarithmic plot of the occupation probability /^ gt in the condensed 
phase (6 = 4, /9 = 4/3c = 2), against the ratio k/M. Top to bottom: M = 20, 
40, 80, 160, and 320. Full lines: fk,st obtained by (5.4) and (5.5). Full (empty) 
symbols: maxima (minima) of occupation probability. Dashed vertical lines: 
asymptotic locations of the minima: k/M = A/(2M) = (p — pc)/2 = 3/4, and 
of the maxima: k/M = A/M = p — pc = 3/2. 



We first analyse the behaviour of the occupation probability fk^st in the condensed 
phase. Figure 3 shows a logarithmic plot of fk^st, computed using equations (5.4) 
and (5.5), against the ratio k/M, for 6 = 4, p = 2, and several values of M. This 
plot exhibits the following features. 

• For k <C M, the distribution fk^st is approximately given by the power law (6.4) of 
an infinite critical system. 
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• The contribution of the condensate slowly builds up as a probability hump around 
k = A = M{p — pc), the mean number of excess particles. 

• One observes a broad and shallow probability "dip" in the region located between 
the critical background and the condensate hump, i.e., in the region A; ^ 1 and 
A-fc> 1. 

The region of the dip is dominated by configurations where the excess particles are 
shared by two sites. Indeed, one has (see [19] for a proof) 

/fc,,t ^{b- i)rW p(/_^), (fc » 1, A - /. » 1). (9.1) 

The observed locations of the maxima {k ^ A) and minima {k ^ A/2) of the occupation 
probability corroborate this picture, as explained in the caption of figure 3. 

These observations lead to the following crude estimate for the characteristic time: 

r ~ -^, (9.2) 

/min 

since the minimum /mm of fk^st corresponds to a barrier to cross, in the spirit of the 
Arrhenius law. The limiting scale of time is that required for the passage of this 
potential barrier. Eq. (9.1) implies that /min is reached near the middle of the dip 
region {k ^ N/2), and therefore (9.2) yields 

T ~ A^ (9.3) 

We now present a more precise treatment. Assume that the condensate is on site 
number 1 at the initial observation time (t = 0). The number A'^i(t) of particles on that 
site is initially very large, A'^i(O) ~ A, and therefore evolves slowly, until the condensate 
dissolves into the critical background. Thus 

• We single out A''i(t) as the collective co-ordinate of the system, that is the 
appropriate slow variable describing the dynamics of the condensate. 

• We model the dynamics of Ni{t) by (4.1), i.e., by a biased diffusive motion on the 
interval k = 0, . . . ,N. The left hopping rate is taken equal to the microscopic rate: 
/ifc = Uk- The right hopping rate Afc is chosen such that, in the stationary state, the 
probability fk^st of the effective model coincide with the occupation probability (5.5) 
of the original ZRP. The detailed balance condition (4.2) yields 

Afc = = ^ = Mst(M - 1, A^ - A;), 

where the right side of the equation is the average rate coming from M — 1 sites, 
containing N — k particles (see (5.11)). The rate A^ is thus a function of k, M, and 
A^. For M = 2, this formula gives Xk = UN-k, as expected. In the fluid phase, in 
the thermodynamic limit, the rates A^ converge to Zq, defined in (5.8). Finally, for 
the condensed phase, in the dip region, we obtain 

Afc ^ 1 + . _ , = MA-fc- 
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This effective description reduces the full model to a Markovian model for one degree 
of freedom in an asymmetric potential. The two valleys of the potential are separated 
by a high (power-law) barrier. The left potential valley, corresponding to the critical 
background, has a weight Pl ~ 1? whereas the right potential valley, corresponding to 
the hump of the condensate, has a weight Pr ~ 1/M <C 1 (see section 6.2). 

In this framework the stationary dynamics of the condensate is characterised by 
a single diverging timescale. We choose to define this timescale, denoted by TMarkov, to 
be the crossing time Tl from the right valley to the left one in the effective Markovian 
problem. The characteristic time is thus expressed by 

Af 1 AT 

■^Markov = -^ L ^ / j 7 / j Jm,sty W""^/ 

in terms of known quantities, the rates fik and the stationary probabilities /fc,st- Its 

asymptotic growth is easily determined by noting that (9.4) is dominated by the 

behaviour of the probability fk^st in the region of the dip. Hence, inserting the 
expression (9.1) into (9.4), and evaluating the sum as an integral, we obtain 

br{b + i) A'^^' _ bTib + i) ,^ ,,^, . 



{b-l)T{2b + 2) M (6-l)r(26 + 2)^^ ^^^ ^^ ' ^^'^^ 

In order to compare the above theoretical predictions to the measured flipping time 
r, we compute the two-time stationary correlation function 

C(t,0) = (iVi(t)iVi(0))-p2. 

This quantity decays exponentially with a relaxation constant which gives a natural 
measure of r. It is found that r ~ TMarkov ~ M'^ in the MF and IDAS geometries, and 
that r ~ MrMarkov ~ M''"'"^ in the IDS geometry. For the latter case the occurence of 
one supplementary power in the system size has the same origin as for nonequilibrium 
dynamics. 

9.4- Last remarks 

Table 1 summarises the values of the dynamical exponents z and Z, such that Tnon-st ~ 
M^ and r ~ M^, where r is the characteristic timescale for the stationary motion of the 
condensate. As recalled above, the non- stationary dynamical exponents are insensitive 



Geometry 


z 


Z 


MF, IDAS 


2 


b 


IDS 


3 


6+1 



Table 1. Non-stationary and stationary dynamical exponents of the ZRP with 
static exponent b > 2. 



to the exponent b, and more generally to the statics, provided the system is in its 
condensed phase. This feature is easily understood in the context of the Markovian 
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Ansatz proposed in the present work. Indeed the last stage of the formation of the 
condensate, i.e., the disappearance of the smaller of the last two precursors, implies no 
barrier crossing. In terms of the occupation of the condensate, it corresponds to the 
transition from A^^i to A, where the initial occupation A^i of the larger precursor was 
already larger than A/2, corresponding to the top of the potential barrier. This explains 
why Tnon-st IS giveu by the diffusive timescale, both in the framework of the Markovian 
Ansatz and in the MF and IDAS geometries. 

10. Further references 

Complementary aspects to the present notes can be found in [11, 38, 39, 40]. 
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Appendix: Proof of eqs. (1.12) and (1.13) 

We recall the fundamental equation (1.16) 

p{Wk,i - Wk,o) + q{Wi,, - Wi,o) = 

P (w,,,,., P-^±^^ - W,o] + g fw^,.M-i ^^^^^ - W,,o] . (A.l) 

V PkPi J \ PkPi J 

We first prove (1.12) for the symmetric case, p = 1/2. Setting Xk = Wk^iPkPi, where 
k + I = n, (A.l) can be rewritten as 

•^k •^n~{k—l) ■^k+l •^n—k- 

This expression is therefore a constant independent of fc, which is equal to zero, as can 
be seen by taking k = n. We thus obtain 

•^fc+l ■^n—k 

which is the detailed balance condition 

PkPi Wk,i = Pk-iPi+i Wi+i^k-i- (A. 2) 

We now show that in the general case, p ^ 1/2, eq. (A.l) yields two constraints on 
the rate: eq. (A. 2) to be interpreted as the pairwise balance condition, and eq. (1.13) 

Wi^k - Wk^i = Wi^^ - Wk,^. (A.3) 

Set ttk = PkPi Wk^i- Eq. (A.l) can be rewritten as 

Vk+i - Vk = {p - q){an~k - ak) (A.4) 
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where 

yk=pxk-qxn-{k-i), 1/0 = 0. (A. 5) 

If 

yk+i = yn^k (A. 6) 

then it follows immediately that Xk+i = Xn-k, which is the condition for pairwise balance 
seen above. This relation itself plugged into (A. 5), yields 

Xk+l ~ Xk = CLn-k ~ O^k 

which is (A. 3). In order to prove (A. 6) we set 

Ak = ai + ... + ak 
i.e. ak = Ak — Ak-i. We thus have 

Vk + {P- q){Ak-l + An-k) = Vk+l + {P- q){Ak + An-k-l) 

which is equal to {p — q')A„, hence 

Vk = {p- q){An - An-k - Ak + ak). 
Therefore yk — {p — q)0'k is symmetric in the change k —>■ n — k, and finally 

Vk- {p- q)ak = yn^k - {p- q)an-k 

which, using (A. 4), yields (A. 6). 
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